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WD + MS systems as the progenitor of SNe la 
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Abstract The single-degenerate (SD) channel for the 
progenitors of type la supernovae (SNe la) is one of 
the most popular channels, in which a carbon-oxygen 
white dwarf (CO WD) accretes hydrogen-rich material 
from its companion, increases its mass to the Chan- 
drasekhar mass limit, and then explodes as a SN la. 
We show the initial and final parameter space for SNe 
la in a (logP 1 , M\) plane and find that the positions of 
some famous recurrent novae, as well as a supersoft X- 
ray source (SSS), RX J0513. 9-6951, are well explained 
by our model. The model can also explain the space 
velocity and mass of Tycho G, which is now suggested 
to be the companion star of Tycho's supernova. Our 
study indicates that the SSS, V Sge, might be the po- 
tential progenitor of supernovae like SN 2002ic if the 
delayed dynamical-instability model due to Han & Pod- 
siadlowski (2006) is appropriate. Following the work of 
Meng, Chen & Han (2009), wc found that the SD model 
(WD + MS) with an optically thick wind can explain 
the birth rate of supernovae like SN 2006X and repro- 
duce the distribution of the color excess of SNe la. The 
model also predicts that at least 75% of all SNe la may 
show a polarization signal in their spectra. 

Keywords stars: binaries: general — stars: super- 
novae: general— individual(SN 2002ic, SN 2006X)— 
stars: white dwarfs 



1 Introduction 

Type la supernovae (SNe la) play an important role as 
cosmological distance indicators to measure cosmologi- 
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cal parameters (e.g. and A; Riess et al. 1998; Perl- 
mutter ct al. 1999), which have led to the discovery 
of the accelerating expansion of the universe. However, 
the exact nature of SNe la progenitors has not been 
clear (see the reviews by Hillebrandt & Nicmcycr 2000; 
Leibundgut 2000). It is widely accepted that a SN la 
originates from the thermonuclear runaway of a carbon- 
oxygen white dwarf (CO WD) in a binary system. The 
CO WD accretes material from its companion, increases 
mass to its maximum stable mass, and then explodes 
as a thermonuclear runaway (Branch 2004). Among 
all the progenitor models, the single-degenerate (SD) 
Chandrasekhar model is widely studied (Yungelson et 
al. 1995; Li & van den Heuvel 1997; Hachisu et al. 
1999a, 1999b; Nomoto ct al. 1999; Langer et al. 2000; 
Han & Podsiadlowski 2004; Chen & Li 2007; Han 2008; 
Meng, Chen & Han 2009; Lii ct al. 2009) and upheld 
by many observation, e.g. the discovery of the potential 
companion of Tycho's supernova (Ruiz-Lapucnte et al. 
2004) and the variable circumstcllar absorption lines 
(Patat et al. 2007). In the SD model, the maximum 
stable mass of the CO WD is ~ 1.378M Q (close to the 
Chandrasekhar mass, Nomoto, Thielemann & Yokoi 
1984), and the companion is probably a main sequence 
star or a slightly evolved star (WD+MS), or a red-giant 
star (WD+RG) (Whclan & Ibcn 1973; Nomoto, Thiele- 
mann & Yokoi 1984). In the paper, we focus on the 
WD+MS channel, which is a very important channel 
to contribute to SNe la. 

The paper is organized as follows: we describe our 
methods in Section 2, show our results in Section 3, 
present some discussions in Section 4, and then finally 
in Section 5, we give our conclusions. 



2 Methods 

Our method for treating the binary evolution of WD 
+ MS systems is the same as that in Meng, Chen & 
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Han (2009), and here, we only give a simply descrip- 
tion. We use the stellar evolution code of Egglcton 
(1971, 1972, 1973) to calculate the binary evolutions of 
WD+MS systems. Instead of solving the stellar struc- 
ture equations of a WD, we adopt the prescription of 
Hachisu et al. (1999a) on WDs accreting hydrogen-rich 
material from their companions. We assume that if the 
mass-transfer rate, |M 2 |, between WD and its compan- 
ion is larger than a critical value, Af cr , the accreted 
hydrogen-rich material steadily burns on the surface of 
WD at the rate of M CI , while the unprocessed matter 
is lost from the system as an optically thick wind at 
a rate of M w i n d = \M%\ — M cr (Hachisu et al. 1996). 
The material lost as the optically thick wind may ex- 
ist as circumstcllar material (CSM) and be a possible 
origin of color excess of SNe la (Meng et al. 2009). 
When M2I locates in the range of Af cr to |Af cr , the 
hydrogen burning is stable or weakly unstable and no 
material is lost from the system. If M2I is smaller 
than ^Af cr , the hydrogen burning is heavily unstable 
and no material accumulates on the surface of the WD. 
The hydrogen-rich material accumulated is converted 
into helium. When the mass of the helium reaches a 
certain value, a helium flash may occur and a part of 
helium may be lost from the system. The mass accu- 
mulation efficiency for helium-shell flashes is from Kato 
& Hachisu (2004). If the mass of a WD can reach 1.378 
M©, we assume that the WD explodes as a SN la. 

3 Results 

3.1 initial and final parameters space 

We carried out a detailed binary evolution calculation 
for more than 25,000 WD+MS systems with various 
mctallicities, and we obtained a large, dense model 
grid. The final outcomes of all the binary evolution cal- 
culations are summarized in the initial orbital period- 
secondary mass (logP 1 , AP,) plane (see Figs. 2 and 3 in 
Meng, Chen & Han 2009). 

Our results may provide help to judge whether a WD 
+ MS system may explode as a SN la. In Fig. [U we 
show the initial contour for SNe la and the final state 
of binary evolution in the (log P — A/2) plane at the mo- 
ment of SNe la explosion for the case of M$y D = 1.2M & 
with Z = 0.02. The cases of other WD masses and 
other mctallicities arc similar. In the figure, we also 
show the position of a supersoft X-ray source (SSS, RX 
J0513. 9-6951) and three famous recurrent novae, i.e. 
U Sco, V394 CrA and CI Aql. Considering the WD 
masses of RX J0513. 9-6951 and U Sco and their po- 
sitions in the (logP — A/2) plane, they are very likely 
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Fig. 1 Parameter regions producing SNe la in the (log P — 
M2) (orbital period-donor mass) plane for the WD + MS 
systems. The initial WD mass is 1.2OM . The WD + MS 
system insides the region encircled by solid line (labeled 
"initial") will increase its white dwarf mass up to 1.378 Mq, 
where we assume a SN la explosion. The final state of the 
WD + MS system in the plane is encircled by a dot-dashed 
line (labeled "final"). Filled squares indicate SN la explo- 
sions during an optically thick wind phase. Filled circles and 
filled triangles denote SN la explosions after the wind phase 
while hydrogen-shell burning is stable or mildly unstable, 
respectively. A supersoft X-ray source, RX J0513. 9-6951, 
(open star) is plotted. Three recurrent novae are indicated 
by a filled star, a solar symbol and an earth symbol, respec- 
tively. 

candidates for the progenitors of SNe la. Due to the un- 
certainties of the companions of V394 CrA and CI Aql, 
they are potential candidates for the progenitors of SNe 
la (See Hachisu et al. (2008) and Meng, Yang & Ceng 
(2009) for similar discussions about these objects). 

3.2 companion properties 

We also show the parameter spaces for companions at 
the moment of supernova explosion. In Fig. [2j we show 
the final states of the companions with various mctal- 
licities at the moment of supernova explosion in the 
(A/f N — V OI b) (final companion mass - orbital velocity) 
plane. We noticed that the suggested companion star of 
Tycho's supernova, Tycho G, locates in the permitted 
region and all mctallicities seem to account for the posi- 
tion of the potential companion. Actually, the detailed 
binary population synthesis results show that only the 
case of Z = 0.02 is appropriate for Tycho G (see also 
Meng, Yang & Geng 2009). 
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Fig. 2 The final states of companion stars at the moment 
of SNe la explosion for different metallicities. The initial 
WD mass is 1.20Mq. Cross represents Tycho G, which is a 
potential candidate for the companion of Tycho's supernova 
(Ruiz-Lapuente et al. 2004; Branch 2004). The length of 
the cross represents observational errors. 



3.3 SN 2002ic 

SN 2002ic is the first object showing hydrogen lines in 
its spectrum, which is interpreted as the interaction be- 
tween supernova ejecta and CSM (Hamuy et al. 2003). 
Many models have been suggested to explain this rare 
object. Among all the models, the delayed dynamical- 
instability model suggested by Han & Podsiadlowski 
(2006) is more interesting and some predictions from 
the model seem to be consistent with observations, es- 
pecially in the sense of the birth rate and delay time 
of the rare object (Aldcring et al. 2006; Prieto et al. 
2007). In the scenario of Han & Podsiadlowski (2006), 
SN 2002ic may be from the WD + MS channel, where 
the CO WD accretes material from its relatively mas- 
sive companion (~ 3.0-Mq), and increases its mass to 
1.30 Mq before experiencing a delayed dynamical in- 
stability. Following the study of Han & Podsiadlowski 
(2006), Meng, Chen & Han (2009) showed the param- 
eter space for SN 2002ic with various metallicities and 
predicated that supernovae like SN 2002ic may not be 
found in extremely low-metallicity environments. How- 
ever, more evidence is necessary to confirm the scenario 
suggested by Han & Podsiadlowski (2006), especially to 
find a progenitor system. 

V Sgc is a well observed quasi-periodic transient SSS 
in our Galaxy. We found that it just locates at the 
boundary of parameters for SN 2002ic in (log P'jMj) 
plane. The mass-loss rate for V Sge can be as large as 
~ 10 MqYI , indicated by radio observation (Lock- 
ley et al. 1997, 1999), which is consistent with the pre- 
diction from the delayed dynamical-instability model. 



So, it is possible that V Sge is the first candidate of 
the progenitor of supernovae like SN 2002ic if the de- 
layed dynamical-instability model in Han & Podsiad- 
lowski (2006) is appropriate. Based on the prediction 
from Han & Podsiadlowski (2006) that 1 in 100 SNe 
la belongs to the subgroup of 2002ic-like supernovae 
and considering that the life time of V Sge is ~ 10 5 
yr, there should be several V Sge-type stars belong- 
ing to the WD + MS system in our Galaxy if we take 
the Galactic birth rate of SNe la as 3-4xl0 _3 yr _1 (van 
den Bergh & Tammann 1991; Cappellaro & Turatto 
1997). Observationally, Steiner & Diaz (1998) listed 
four V Sge-type stars in the Galaxy and discussed their 
similar spectroscopic and photometric properties. So, 
the number of V Sge-type stars is also consistent with 
the prediction from the delayed dynamical-instability 
model (see Meng, Yang & Gcng 2009 for details). 

However, in our model grid, the case of M^-q = 
\.2Mq is the only one to account for the position of V 
Sge in the (logP — M2) plane, but it currently seems 
that 1.2Mq CO WDs are probably only in special cir- 
cumstance to form such as an extremely high metal- 
licity (Meng, Chen & Han 2008). In addition, the CO 
WD with a mass high than 1.1M seems not to explode 
as a SN la (Umeda et al. 1999). So it is premature to 
obtain a definitive conclusion that there is a relation be- 
tween the V Sge-type stars and the delayed dynamical 
instability model. 

3.4 the birth rate 

Incorporating our results into Hurleys rapid stellar evo- 
lution code (Hurley et al. 2000, 2002), we calculate the 
birth rate of SNe la for the WD+MS model. We use 
Monte Carlo simulation to generate a primordial binary 
sample (see Meng, Chen & Han 2009 for the basic pa- 
rameters for Monte Carlo simulation). We found that 
for the case of single star burst, most supernovae occur 
between 0.1 Gyr and 2 Gyr after star formation, and a 
high mctallicity leads to a systematically earlier explo- 
sion time. The peak value of the birth rate increases 
with metallicity Z . However, the WD+MS model only 
can account for about 30 percent of the Galactic birth 
rate (see also Meng, Chen & Han 2009). 

3.5 SN 2006X 

SN 2006X is the first case to show the signal of cir- 
cumstellar material, which shows its single-degenerate 
nature (Patat et al. 2007). Patat et al. (2007) sug- 
gested that the progenitor of SN 2006X is a WD + RG 
system based on the expansion velocity of the circum- 
stellar material. Blondin et al. (2009) found that 6 in 
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100 SNe la should belong to the rare group. Hachisu 
ct al. (2008) suggested that SN 2006X may be from a 
WD + MS system, and if a WD explodes at the opti- 
cally thick wind phase, the SN may show a signal like 
SN 2006X. Following the study of Meng, Chen & Han 
(2009), we calculated the evolution of the birth rate of 
supernovae like SN 2006X by assuming that if a SN 
la explodes at the optically thick wind phase, it is a 
2006X-like supernova. We found that the progenitor 
age of SN 2006X is between 0.3Gyr and lGyr, which 
means that there is star formation during the recent 
1 Gyr in the host galaxy of supernovae like SN 2006X. 
The host galaxies of SN 2006X and its twins, SN 1999cl, 
are both spiral galaxies, and show a significant signal 
of star formation at present (Kanpen et al. 1993, 1996; 
Wong & Blitz 2002). We also noticed that supernovae 
like SN 2006X are a relatively rare subclass of SNe la: 1 
to 14 in 100 SNe la can be of this type, which depends 
on the common envelope ejection efficiency, acE- The 
birth rate found by Blondin et al. (2009) if located 
in the range. So, our results uphold the suggestion of 
Hachisu ct al. (2008) that the progenitor system of a 
supernova like SN 2006X may be a WD + MS system 
(see Meng, Yang & Ceng 2009b). 

3.6 color excess 

In our binary evolution calculation, we assume that a 
part of the hydrogen-rich material may be lost from 
the binary system as an optically thick wind. The lost 
material may exist as circumstcllar material (CSM). If 
this scenario is appropriate, the CSM should be the ori- 
gin of the color excess of SNe la. Reindl ct al. (2005) 
showed the color excesses to be more than one hundred 
SNe la at maximum light, which provides an opportu- 
nity to check whether the SD model may reproduce the 
distribution of the color excess of SNe la. Following the 
study of Meng, Chen & Han (2009) and via a simple 
analytic method, we may reproduce the distribution of 
color excesses of SNe la by our binary population syn- 
thesis (BPS) approach if the velocity of the optically 
thick wind is taken to be of the order of 10 km s _1 . 
However, if the wind velocity is larger than 100 km 
s _1 , the reproduction is bad (see Meng et al. 2009 for 
details). 

3.7 Tycho G 

Following the study of Meng, Chen & Han (2009) and 
via a BPS approach, we calculated the distributions of 
the parameters of the binary systems at the moment 
of supernova explosion and the properties of compan- 
ions after supernova explosion, e.g. mass, radius, space 



velocity and surface gravity. The former may provide 
physics input when one simulates the interaction be- 
tween supernova ejecta and its companion, and the lat- 
ter may help to search for the companions in supernova 
remnants. It is surprising that the properties of the 
potential companion of Tycho's supernova, Tycho G, 
are well consistent with our BPS results, and only the 
case of Z = 0.02 may account for the properties of Ty- 
cho G (please, see our following paper: Meng & Yang 
2009). However, a conclusive result about Tycho G is 
premature, since there have been considerable debates 
about whether the star is related to Tycho's supernova 
(Kerzendorf et al. 2009; Hernandez ct al. 2009). 

3.8 polarization 

After the explosion, the explosion ejecta impacts into 
the envelope of the companion and strips off some 
hydrogen- rich material from the surface of the compan- 
ion. Because of the existence of the companion, a hole 
forms in the explosion ejecta and the hole will never 
disappear since the movement of the ejecta is super- 
sonic. Then the explosion remnant is asphcric, and the 
aspheric structure may reveal itself by a polarized spec- 
trum (Marietta et al. 2000; Meng, Chen & Han 2007). 
Based on the numerical simulation of Marietta et al. 
(2000) and Kasen et al. (2004), we found that at least 
75% of all SNe la can be detected by observation of 
the polarization. At present, all supernovae which are 
observed by spcctropolarimetry have various degrees of 
polarization signal (Leonard ct al. 2005). 

4 Discussion 

From our study, the WD + MS channel can only ac- 
count for about 1/3 of the SNe la observed. Therefore, 
there may be other channels or mechanisms contribut- 
ing to SNe la. A wide symbiotic system, WD + RG, is 
a possible progenitor of SNe la. Now, we are construct- 
ing a complete model including WD + MS channel and 
WD + RG channel, hoping to resolve the problem. In 
our new model, the range of the companion mass for 
SNe la is from 0.6 M Q to 5 M . The mass of WD 
contributing to SNe la may be as low as 0.565 M Q . 

An alternative is the double-degenerate (DD) chan- 
nel. Although, it is theoretically less favored, observa- 
tion showed that it may contribute to a few SNe la, e.g. 
SN 2003fg and SN 2006gz (Howell ct al. 2006; Hicken 
et al. 2007). So, at present it is premature to obtain a 
definitive conclusion for the DD model. 
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5 Conclusion 



The WD+MS channel is a very important channel for 
SNe la and the channel may contribute to SNe la by 
about 30%. Based on a detailed binary evolution cal- 
culation, we show the initial and final parameter space 
for SNe la in the (log P'jM^) plane. The results may 
provide help to survey potential candidates of SNe la. 
Following the results of Meng, Chen & Han (2009) and 
via a binary population synthesis approach, we found 
that the SD channel may reproduce the birth rate of 
SN 2006X and the distribution of color excess of SNe 
la. The properties of the potential companion of Ty- 
cho's supernova, Tycho G, may also be consistent with 
our BPS results. 
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